Contribution of Microorganisms with the Clade II Nitrous Oxide Reductase to Suppression of Surface Emissions of Nitrous Oxide.

The sources and sinks of nitrous oxide, as control emissions to the atmosphere, are generally poorly constrained for most environmental systems. Initial depth-resolved analysis of nitrous oxide flux from observation wells and the proximal surface within a nitrate contaminated aquifer system revealed high subsurface production but little escape from the surface. To better understand the environmental controls of production and emission at this site, we used a combination of isotopic, geochemical, and molecular analyses to show that chemodenitrification and bacterial denitrification are major sources of nitrous oxide in this subsurface, where low DO, low pH, and high nitrate are correlated with significant nitrous oxide production. Depth-resolved metagenomes showed that consumption of nitrous oxide near the surface was correlated with an enrichment of Clade II nitrous oxide reducers, consistent with a growing appreciation of their importance in controlling release of nitrous oxide to the atmosphere. Our work also provides evidence for the reduction of nitrous oxide at a pH of 4, well below the generally accepted limit of pH 5.

Seasonality and longer-term development generate temporal dynamics in the Populus microbiome.

Temporal variation in community composition is central to our understanding of the assembly and functioning of microbial communities, yet the controls over temporal dynamics for microbiomes of long-lived plants, such as trees, remain unclear. Temporal variation in tree microbiomes could arise primarily from seasonal (i.e., intra-annual) fluctuations in community composition or from longer-term changes across years as host plants age. To test these alternatives, we experimentally isolated temporal variation in plant microbiome composition using a common garden and clonally propagated plants, and we used amplicon sequencing to characterize bacterial/archaeal and fungal communities in the leaf endosphere, root endosphere, and rhizosphere of two Populus spp. over four seasons across two consecutive years. Microbial community composition differed among seasons and years (which accounted for up to 21% of the variation in microbial community composition) and was correlated with seasonal dissimilarity in climatic conditions. However, microbial community dissimilarity was also positively correlated with time, reflecting longer-term compositional shifts as host trees aged. Together, our findings demonstrate that temporal patterns in tree microbiomes arise from both seasonal fluctuations and longer-term changes, which interact to generate unique seasonal patterns each year. In addition to shedding light on two important controls over the assembly of plant microbiomes, our results also suggest future studies of tree microbiomes should account for background temporal dynamics when testing the drivers of spatial patterns in microbial community composition and temporal responses of plant microbiomes to environmental change.IMPORTANCEMicrobiomes are integral to the health of host plants, but we have a limited understanding of the factors that control how the composition of plant microbiomes changes over time. Especially little is known about the microbiome of long-lived trees, relative to annual and non-woody plants. We tested how tree microbiomes changed between seasons and years in poplar (genus Populus), which are widespread and ecologically important tree species that also serve as important biofuel feedstocks. We found the composition of bacterial, archaeal, and fungal communities differed among seasons, but these seasonal differences depended on year. This dependence was driven by longer-term changes in microbial composition as host trees developed across consecutive years. Our findings suggest that temporal variation in tree microbiomes is driven by both seasonal fluctuations and longer-term (i.e., multiyear) development.

Elevated temperature alters microbial communities, but not decomposition rates, during 3 years of in situ peat decomposition.

IMPORTANCE: Microbial community changes in response to climate change drivers have the potential to alter the trajectory of important ecosystem functions. In this paper, we show that while microbial communities in peatland systems responded to manipulations of temperature and CO2 concentrations, these changes were not associated with similar responses in peat decomposition rates over 3 years. It is unclear however from our current studies whether this functional resiliency over 3 years will continue over the longer time scales relevant to peatland ecosystem functions.

Variation in carbon and nitrogen concentrations among peatland categories at the global scale.

Peatlands account for 15 to 30% of the world's soil carbon (C) stock and are important controls over global nitrogen (N) cycles. However, C and N concentrations are known to vary among peatlands contributing to the uncertainty of global C inventories, but there are few global studies that relate peatland classification to peat chemistry. We analyzed 436 peat cores sampled in 24 countries across six continents and measured C, N, and organic matter (OM) content at three depths down to 70 cm. Sites were distinguished between northern (387) and tropical (49) peatlands and assigned to one of six distinct broadly recognized peatland categories that vary primarily along a pH gradient. Peat C and N concentrations, OM content, and C:N ratios differed significantly among peatland categories, but few differences in chemistry with depth were found within each category. Across all peatlands C and N concentrations in the 10-20 cm layer, were 440 ± 85.1 g kg-1 and 13.9 ± 7.4 g kg-1, with an average C:N ratio of 30.1 ± 20.8. Among peatland categories, median C concentrations were highest in bogs, poor fens and tropical swamps (446-532 g kg-1) and lowest in intermediate and extremely rich fens (375-414 g kg-1). The C:OM ratio in peat was similar across most peatland categories, except in deeper samples from ombrotrophic tropical peat swamps that were higher than other peatlands categories. Peat N concentrations and C:N ratios varied approximately two-fold among peatland categories and N concentrations tended to be higher (and C:N lower) in intermediate fens compared with other peatland types. This study reports on a unique data set and demonstrates that differences in peat C and OM concentrations among broadly classified peatland categories are predictable, which can aid future studies that use land cover assessments to refine global peatland C and N stocks.

Relationships between Sphaerulina musiva Infection and the Populus Microbiome and Metabolome.

Pathogenic fungal infections in plants may, in some cases, lead to downstream systematic impacts on the plant metabolome and microbiome that may either alleviate or exacerbate the effects of the fungal pathogen. While Sphaerulina musiva is a well-characterized fungal pathogen which infects Populus tree species, an important wood fiber and biofuel feedstock, little is known about its systematic effects on the metabolome and microbiome of Populus. Here, we investigated the metabolome of Populus trichocarpa and Populus deltoides leaves and roots and the microbiome of the leaf and root endospheres, phylloplane, and rhizosphere to understand the systematic impacts of S. musiva abundance and infection on Populus species in a common garden field setting. We found that S. musiva is indeed present in both P. deltoides and P. trichocarpa, but S. musiva abundance was not statistically related to stem canker onset. We also found that the leaf and root metabolomes significantly differ between the two Populus species and that certain leaf metabolites, particularly the phenolic glycosides salirepin and salireposide, are diminished in canker-infected P. trichocarpa trees compared to their uninfected counterparts. Furthermore, we found significant associations between the metabolome, S. musiva abundance, and microbiome composition and α-diversity, particularly in P. trichocarpa leaves. Our results show that S. musiva colonizes both resistant and susceptible hosts and that the effects of S. musiva on susceptible trees are not confined to the site of canker infection. IMPORTANCE Poplar (Populus spp.) trees are ecologically and economically important trees throughout North America. However, many western North American poplar plantations are at risk due to the introduction of the nonnative fungal pathogen Sphaerulina musiva, which causes leaf spot and cankers, limiting their production. To better understand the interactions among the pathogen S. musiva, the poplar metabolome, and the poplar microbiome, we collected leaf, root, and rhizosphere samples from poplar trees consisting of 10 genotypes and two species with differential resistance to S. musiva in a common garden experiment. Here, we outline the nuanced relationships between the poplar metabolome, microbiome, and S. musiva, showing that S. musiva may affect poplar trees in tissues distal to the site of infection (i.e., stem). Our research contributes to improving the fundamental understanding of S. musiva and Populus sp. ecology and the utility of a holobiont approach in understanding plant disease.

Functional Redundancy in Soil Microbial Community Based on Metagenomics Across the Globe.

Understanding the contribution of soil microbial communities to ecosystem processes is critical for predicting terrestrial ecosystem feedbacks under changing climate. Our current understanding lacks a consistent strategy to formulate the linkage between microbial systems and ecosystem processes due to the presumption of functional redundancy in soil microbes. Here we present a global soil microbial metagenomic analysis to generalize patterns of microbial taxonomic compositions and functional potentials across climate and geochemical gradient. Our analyses show that soil microbial taxonomic composition varies widely in response to climate and soil physicochemical gradients, while microbial functional attributes based on metagenomic gene abundances are redundant. Among 17 climate zones, microbial taxonomic compositions were more distinct than functional potentials, as climate and edaphic properties showed more significant influence on microbial taxonomic compositions than on functional potentials. Microbial taxonomies formed a larger and more complex co-occurrence network with more module structures than functional potentials. Functional network was strongly inter-connected among different categories, whereas taxonomic network was more positively interactive in the same taxonomic groups. This study provides strong evidence to support the hypothesis of functional redundancy in soil microbes, as microbial taxonomic compositions vary to a larger extent than functional potentials based on metagenomic gene abundances in terrestrial ecosystems across the globe.

Beyond the usual suspects: methanogenic communities in eastern North American peatlands are also influenced by nickel and copper concentrations.

Peatlands both accumulate carbon and release methane, but their broad range in environmental conditions means that the diversity of microorganisms responsible for carbon cycling is still uncertain. Here, we describe a community analysis of methanogenic archaea responsible for methane production in 17 peatlands from 36 to 53 N latitude across the eastern half of North America, including three metal-contaminated sites. Methanogenic community structure was analysed through Illumina amplicon sequencing of the mcrA gene. Whether metal-contaminated sites were included or not, metal concentrations in peat were a primary driver of methanogenic community composition, particularly nickel, a trace element required in the F430 cofactor in methyl-coenzyme M reductase that is also toxic at high concentrations. Copper was also a strong predictor, likely due to inhibition at toxic levels and/or to cooccurrence with nickel, since copper enzymes are not known to be present in anaerobic archaea. The methanogenic groups Methanocellales and Methanosarcinales were prevalent in peatlands with low nickel concentrations, while Methanomicrobiales and Methanomassiliicoccales were abundant in peatlands with higher nickel concentrations. Results suggest that peat-associated trace metals are predictors of methanogenic communities in peatlands.

A k-mer based approach for classifying viruses without taxonomy identifies viral associations in human autism and plant microbiomes.

Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an essential role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus' genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples, even if they do not have a taxonomy. To identify and quantify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses to create a database of metagenomic viruses. The resulting Kraken2 database of the metagenomic viruses can be found here: https://www.osti.gov/biblio/1615774 and is compatible with Kraken2. We then integrated the viral classification database with databases created with genomes from NCBI for use with ParaKraken (a parallelized version of Kraken provided in Supplemental Zip 1), a metagenomic/transcriptomic classifier. To illustrate the breadth of our utility for classifying metagenome viruses, we analyzed data from a plant metagenome study identifying genotypic and compartment specific differences between two Populus genotypes in three different compartments. We also identified a significant increase in abundance of eight viral sequences in post mortem brains in a human metatranscriptome study comparing Autism Spectrum Disorder patients and controls. We also show the potential accuracy for classifying viruses by utilizing both the JGI and NCBI viral databases to identify the uniqueness of viral sequences. Finally, we validate the accuracy of viral classification with NCBI databases containing viruses with taxonomy to identify pathogenic viruses in known COVID-19 and cassava brown streak virus infection samples. Our method represents the compulsory first step in better understanding the role of viruses in the microbiome by allowing for a more complete identification of sequences without taxonomy. Better classification of viruses will improve identifying associations between viruses and their hosts as well as viruses and other microbiome members. Despite the lack of taxonomy, this database of metagenomic viruses can be used with any tool that utilizes a taxonomy, such as Kraken, for accurate classification of viruses.

Heterospecific Neighbor Plants Impact Root Microbiome Diversity and Molecular Function of Root Fungi.

Within the forest community, competition and facilitation between adjacent-growing conspecific and heterospecific plants are mediated by interactions involving common mycorrhizal networks. The ability of plants to alter their neighbor's microbiome is well documented, but the molecular biology of plant-fungal interactions during competition and facilitation has not been previously examined. We used a common soil-plant bioassay experiment to study molecular plant-microbial interactions among rhizosphere communities associated with Pinus taeda (native host) and Populus trichocarpa (non-native host). Gene expression of interacting fungal and bacterial rhizosphere communities was compared among three plant-pairs: Populus growing with Populus, Populus with Pinus, and Pinus with Pinus. Our results demonstrate that heterospecific plant partners affect the assembly of root microbiomes, including the changes in the structure of host specific community. Comparative metatranscriptomics reveals that several species of ectomycorrhizal fungi (EMF) and saprotrophic fungi exhibit different patterns of functional and regulatory gene expression with these two plant hosts. Heterospecific plants affect the transcriptional expression pattern of EMF host-specialists (e.g., Pinus-associated Suillus spp.) on both plant species, mainly including the genes involved in the transportation of amino acids, carbohydrates, and inorganic ions. Alteration of root microbiome by neighboring plants may help regulate basic plant physiological processes via modulation of molecular functions in the root microbiome.

Minnesota peat viromes reveal terrestrial and aquatic niche partitioning for local and global viral populations.

BACKGROUND: Peatlands are expected to experience sustained yet fluctuating higher temperatures due to climate change, leading to increased microbial activity and greenhouse gas emissions. Despite mounting evidence for viral contributions to these processes in peatlands underlain with permafrost, little is known about viruses in other peatlands. More generally, soil viral biogeography and its potential drivers are poorly understood at both local and global scales. Here, 87 metagenomes and five viral size-fraction metagenomes (viromes) from a boreal peatland in northern Minnesota (the SPRUCE whole-ecosystem warming experiment and surrounding bog) were analyzed for dsDNA viral community ecological patterns, and the recovered viral populations (vOTUs) were compared with our curated PIGEON database of 266,125 vOTUs from diverse ecosystems. RESULTS: Within the SPRUCE experiment, viral community composition was significantly correlated with peat depth, water content, and carbon chemistry, including CH4 and CO2 concentrations, but not with temperature during the first 2 years of warming treatments. Peat vOTUs with aquatic-like signatures (shared predicted protein content with marine and/or freshwater vOTUs) were significantly enriched in more waterlogged surface peat depths. Predicted host ranges for SPRUCE vOTUs were relatively narrow, generally within a single bacterial genus. Of the 4326 SPRUCE vOTUs, 164 were previously detected in other soils, mostly peatlands. None of the previously identified 202,371 marine and freshwater vOTUs in our PIGEON database were detected in SPRUCE peat, but 0.4% of 80,714 viral clusters (VCs, grouped by predicted protein content) were shared between soil and aquatic environments. On a per-sample basis, vOTU recovery was 32 times higher from viromes compared with total metagenomes. CONCLUSIONS: Results suggest strong viral "species" boundaries between terrestrial and aquatic ecosystems and to some extent between peat and other soils, with differences less pronounced at higher taxonomic levels. The significant enrichment of aquatic-like vOTUs in more waterlogged peat suggests that viruses may also exhibit niche partitioning on more local scales. These patterns are presumably driven in part by host ecology, consistent with the predicted narrow host ranges. Although more samples and increased sequencing depth improved vOTU recovery from total metagenomes, the substantially higher per-sample vOTU recovery after viral particle enrichment highlights the utility of soil viromics. Video abstract The importance of Minnesota peat viromes in revealing terrestrial and aquatic niche partitioning for viral populations.

Climate-driven divergence in plant-microbiome interactions generates range-wide variation in bud break phenology.

Soil microbiomes are rapidly becoming known as an important driver of plant phenotypic variation and may mediate plant responses to environmental factors. However, integrating spatial scales relevant to climate change with plant intraspecific genetic variation and soil microbial ecology is difficult, making studies of broad inference rare. Here we hypothesize and show: 1) the degree to which tree genotypes condition their soil microbiomes varies by population across the geographic distribution of a widespread riparian tree, Populus angustifolia; 2) geographic dissimilarity in soil microbiomes among populations is influenced by both abiotic and biotic environmental variation; and 3) soil microbiomes that vary in response to abiotic and biotic factors can change plant foliar phenology. We show soil microbiomes respond to intraspecific variation at the tree genotype and population level, and geographic variation in soil characteristics and climate. Using a fully reciprocal plant population by soil location feedback experiment, we identified a climate-based soil microbiome effect that advanced and delayed bud break phenology by approximately 10 days. These results demonstrate a landscape-level feedback between tree populations and associated soil microbial communities and suggest soil microbes may play important roles in mediating and buffering bud break phenology with climate warming, with whole ecosystem implications.

Cultivating the Bacterial Microbiota of Populus Roots.

The integral role of microbial communities in plant growth and health is now widely recognized, and, increasingly, the constituents of the microbiome are being defined. While phylogenetic surveys have revealed the taxa present in a microbiome and show that this composition can depend on, and respond to, environmental perturbations, the challenge shifts to determining why particular microbes are selected and how they collectively function in concert with their host. In this study, we targeted the isolation of representative bacterial strains from environmental samples of Populus roots using a direct plating approach and compared them to amplicon-based sequencing analysis of root samples. The resulting culture collection contains 3,211 unique isolates representing 10 classes, 18 orders, 45 families, and 120 genera from 6 phyla, based on 16S rRNA gene sequence analysis. The collection accounts for ∼50% of the natural community of plant-associated bacteria as determined by phylogenetic analysis. Additionally, a representative set of 553 had their genomes sequenced to facilitate functional analyses. The top sequence variants in the amplicon data, identified as Pseudomonas, had multiple representatives within the culture collection. We then explore a simplified microbiome, comprised of 10 strains representing abundant taxa from environmental samples, and tested for their ability to reproducibly colonize Populus root tissue. The 10-member simplified community was able to reproducibly colonize on Populus roots after 21 days, with some taxa found in surface-sterilized aboveground tissue. This study presents a comprehensive collection of bacteria isolated from Populus for use in exploring microbial function and community inoculation experiments to understand basic concepts of plant and environmental selection. IMPORTANCE Microbial communities play an integral role in the health and survival of their plant hosts. Many studies have identified key members in these communities and led to the use of synthetic communities for elucidating their function; however, these studies are limited by the available cultured bacterial representatives. Here, we present a bacterial culture collection comprising 3,211 isolates that is representative of the root community of Populus. We then demonstrate the ability to examine underlying microbe-microbe interactions using a synthetic community approach. This culture collection will allow for the greater exploration of the microbial community function through targeted experimentation and manipulation.

Soil metabolome response to whole-ecosystem warming at the Spruce and Peatland Responses under Changing Environments experiment.

In this study, a suite of complementary environmental geochemical analyses, including NMR and gas chromatography-mass spectrometry (GC-MS) analyses of central metabolites, Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) of secondary metabolites, and lipidomics, was used to investigate the influence of organic matter (OM) quality on the heterotrophic microbial mechanisms controlling peatland CO2, CH4, and CO2:CH4 porewater production ratios in response to climate warming. Our investigations leverage the Spruce and Peatland Responses under Changing Environments (SPRUCE) experiment, where air and peat warming were combined in a whole-ecosystem warming treatment. We hypothesized that warming would enhance the production of plant-derived metabolites, resulting in increased labile OM inputs to the surface peat, thereby enhancing microbial activity and greenhouse gas production. Because shallow peat is most susceptible to enhanced warming, increases in labile OM inputs to the surface, in particular, are likely to result in significant changes to CO2 and CH4 dynamics and methanogenic pathways. In support of this hypothesis, significant correlations were observed between metabolites and temperature consistent with increased availability of labile substrates, which may stimulate more rapid turnover of microbial proteins. An increase in the abundance of methanogenic genes in response to the increase in the abundance of labile substrates was accompanied by a shift toward acetoclastic and methylotrophic methanogenesis. Our results suggest that as peatland vegetation trends toward increasing vascular plant cover with warming, we can expect a concomitant shift toward increasingly methanogenic conditions and amplified climate-peatland feedbacks.

Assembly of the Populus Microbiome Is Temporally Dynamic and Determined by Selective and Stochastic Factors.

Recent work shows that the plant microbiome, particularly the initial assembly of this microbiome, influences plant health, survival, and fitness. Here, we characterize the initial assembly of the Populus microbiome across ten genotypes belonging to two poplar species in a common garden using 16S rRNA gene and ITS2 region amplicon sequencing of the leaf endosphere, leaf surface, root endosphere, and rhizosphere. We sampled these microbiomes three times throughout the first growing season and found that the composition of the microbiome changed dramatically over time across all plant-associated habitats and host genotypes. For archaea and bacteria, these changes were dominated by strong homogenizing selection (accounting for 29 to 62% of pairwise comparisons). However, fungal assembly was generally characterized by multiple ecological assembly processes (i.e., a mix of weak selective and dispersal processes). Interestingly, genotype, while a significant moderator of microbiome composition, generally explained less variation than sample date across plant-associated habitats. We defined a set of core genera that accounted for, on average, 36% of the microbiome. The relative abundance of this core community was consistent over time. Additionally, using source tracking modeling, we determined that new microbial taxa colonize from both aboveground and belowground sources, and combined with our ecological assembly null models, we found that both selective and dispersal processes explained the differences between exo- (i.e., leaf surface and rhizosphere) and endospheric microbiomes. Taken together, our results suggest that the initial assembly of the Populus microbiome is time-, genotype-, and habitat-dependent and is moderated by both selective and stochastic factors. IMPORTANCE The initial assembly of the plant microbiome may establish the trajectory of forthcoming microbiome states, which could determine the overall future health of the plant. However, while much is known about the initial microbiome assembly of grasses and agricultural crops, less is known about the initial microbiome of long-lived trees, such as poplar (Populus spp.). Thus, a greater understanding of initial plant microbiome assembly in an ecologically and economically important plant such as Populus is highly desirable. Here, we show that the initial microbiome community composition and assembly in the first growing season of Populus is temporally dynamic and is determined by a combination of both selective and stochastic factors. Our findings could be used to prescribe ecologically informed microbial inoculations and better predict the composition of the Populus microbiome into the future and to better understand its influence on plant health.

Fire alters plant microbiome assembly patterns: integrating the plant and soil microbial response to disturbance.

It is increasingly evident that the plant microbiome is a strong determinant of plant health. While the ability to manipulate the microbiome in plants and ecosystems recovering from disturbance may be useful, our understanding of the plant microbiome in regenerating plant communities is currently limited. Using 16S ribosomal RNA (rRNA) gene and internal transcribed spacer (ITS) region amplicon sequencing, we characterized the leaf, stem, fine root, rhizome, and rhizosphere microbiome of < 1-yr-old aspen saplings and the associated bulk soil after a recent high-intensity prescribed fire across a burn severity gradient. Consistent with previous studies, we found that soil microbiomes are responsive to fire. We extend these findings by showing that certain plant tissue microbiomes also change in response to fire. Differences in soil microbiome compositions could be attributed to soil chemical characteristics, but, generally, plant tissue microbiomes were not related to plant tissue elemental concentrations. Using source tracking modeling, we also show that fire influences the relative dominance of microbial inoculum and the vertical inheritance of the sapling microbiome from the parent tree. Overall, our results demonstrate how fire impacts plant microbiome assembly, diversity, and composition and highlights potential for further research towards increasing plant fitness and ecosystem recovery after fire events.

Phylogenetic diversity of 200+ isolates of the ectomycorrhizal fungus Cenococcum geophilum associated with Populus trichocarpa soils in the Pacific Northwest, USA and comparison to globally distributed representatives.

The ectomycorrhizal fungal symbiont Cenococcum geophilum is of high interest as it is globally distributed, associates with many plant species, and has resistance to multiple environmental stressors. C. geophilum is only known from asexual states but is often considered a cryptic species complex, since extreme phylogenetic divergence is often observed within nearly morphologically identical strains. Alternatively, C. geophilum may represent a highly diverse single species, which would suggest cryptic but frequent recombination. Here we describe a new isolate collection of 229 C. geophilum isolates from soils under Populus trichocarpa at 123 collection sites spanning a ~283 mile north-south transect in Western Washington and Oregon, USA (PNW). To further understanding of the phylogenetic relationships within C. geophilum, we performed maximum likelihood and Bayesian phylogenetic analyses to assess divergence within the PNW isolate collection, as well as a global phylogenetic analysis of 789 isolates with publicly available data from the United States, Japan, and European countries. Phylogenetic analyses of the PNW isolates revealed three distinct phylogenetic groups, with 15 clades that strongly resolved at >80% bootstrap support based on a GAPDH phylogeny and one clade segregating strongly in two principle component analyses. The abundance and representation of PNW isolate clades varied greatly across the North-South range, including a monophyletic group of isolates that spanned nearly the entire gradient at ~250 miles. A direct comparison between the GAPDH and ITS rRNA gene region phylogenies, combined with additional analyses revealed stark incongruence between the ITS and GAPDH gene regions, consistent with intra-species recombination between PNW isolates. In the global isolate collection phylogeny, 34 clades were strongly resolved using Maximum Likelihood and Bayesian approaches (at >80% MLBS and >0.90 BPP respectively), with some clades having intra- and intercontinental distributions. Together these data are highly suggestive of divergence within multiple cryptic species, however additional analyses such as higher resolution genotype-by-sequencing approaches are needed to distinguish potential species boundaries and the mode and tempo of recombination patterns.